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ABSTRACT 

We place constraints on the matter density of the Universe and the amphtude 
of clustering using measurements of the galaxy two-point correlation function from 
the Sloan Digital Sky Survey (SDSS). We generate model predictions for different 
cosmologies by populating rescaled N-body simulations with galaxies using the subhalo 
abundance matching (SHAM) technique. We find = 0.29 ±0.03 and as = 0.86 ± 
0.04 at 68% confidence from fitting the observed two-point galaxy correlation function 
of galaxies brighter than Mr = -18 in a volume limited sample of galaxies obtained by 
the SDSS. We discuss and quantify potential sources of systematic error, and conclude 
that while there is scope for improving its robustness, the technique presented in this 
paper provides a powerful low redshift constraint on the cosmological parameters that 
is complementary to other commonly used methods. 



1 INTRODUCTION 



Determining cosmological parameters such as the matter 
density of the Universe and the amplitude of clustering is 
a key goal of cosmology. Measurements of the clustering of 
galaxies provide a unique window into the distribution of 
dark matter in the Universe from which cosmological pa- 
rameters can be inferred. Galaxy surveys map the distribu- 
tion of visible baryons which indirectly trace the underlying 
distribution of the gravitationally dominant dark matter. 
While the distribution of dark matter for a given cosmology 
can be reliably computed from first principles using cosmo- 
logical N-body simulations, because of uncertainties in the 
detailed physics of galaxy formation - gas cooling, star for- 
mation, feedback etc., the distribution of galaxies cannot be 
robustly predicted from first principles. Consequently, one 
of the principal impediments in inferring cosmological in- 
formation from observations of galaxy clustering is galaxy 
bias, the difference between the distribution of galaxies and 
the underlying dark matter. In this paper, we use a simple 
non-parametric model to relate observed galaxy luminosity 
to halo mass in an N-body simulation, and use this model 
to place constraints on the universal matter density (Qm) 
and the amplitude of clustering (as) by fitting the observed 
clustering of galaxies. 

We relate (sub)halos in our N-body simulation to ob- 
served galaxies using subhalo abundance matching (SIfAM). 
SHAM is a simple, non-parametric model based on assuming 
a monotonic relationship between observed galaxy luminos- 
ity and simulated (sub)halo mass. In this model, all galaxies 
are assumed to be contained within dark matter subhalos, 
and galaxy luminosity is assumed to be monotonically re- 



lated to the present day subhalo mass for central subhalos 
and the subhalo mass at the accretion epoch for satellite 
subhalos. SHAM takes only the space density of galaxies as 
input and predicts the clustering of a galaxy population. We 
use it to predict the clustering of observed volume limited 
samples of galaxies. Distinct cosmologies pro duce distinct 
popu lations of dark matter halos and subhalos jZheng et al.l 
'2002), thus for each cosmology, the SHAM model predicts a 
distinct galaxy auto correlation function. 

SHAM has been successfully used to match theoretica l 
predictions to observables. For example, IGuo et al] (|201G| ) 
match the observed stellar mass - halo mass relation by 
populating N-body simulat ions using the SHAM technique. 
iTruiillo-Gomez et al.l (|201ll 'l match the observed circular ve - 
locity statistics using a similar method. ISimha et al.l (|2012l ) 
found reasonable agreement between its assumptions and 
the output from a cosmological smoothed particle hydrody- 
namics (SPH) simulation which incorporated gas physics, 
star formation and feedback from supernovae driven winds. 

Running a suite of high resolution cosmological N- 
body simulations to adequately sample the cosmological pa- 
rameter space is computationally demanding. Instead, we 
rescale the masses, positions and velocities of halos/subhalos 
obtained from one simulation to different target cosmolo- 
gi es. We use the reduced v ersion of the rescalin g technique 
of lAngulo fc Whit"3 (|2010l 'l as implemented bv iRuiz et al.l 
to construct su bhalo catalogues fo r a giv e n set of cos- 
molog i cal parame t ers. Angulo fc White! l|2010l 'l ; iRuiz et al] 
('201lf) ; iGuo et al.l (|2012l ') have shown that the relevant prop- 
erties of (sub)halos in such scaled models are very close to 
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those in simulations carried out directly with the target cos- 
mology. 

Our approa c h is sim ilar in spirit to previous studies such 
as [Tinker et all l|2012l ) who use halo occupation distribu- 
tion (HOD) models, and lHarker et all (|2007l ) who use semi- 
analytic models and a similar simulation rescaling technique 
to fit the observed clustering of galaxies to infer cosmolog- 
ical parameters. In contra st to the semi - analyt ic models of 
galaxy formation used by iHarker et all Ho^, we do not 
make any assumptions about the detailed gas physics that is 
involved in the process of gala xy formation. Although HOD 
models such as those used bv [Tinker et al.l (|2012l ) are also 
based on statistical descriptions rather than detailed mod- 
elling of the gas physics, there are significant difTerences be- 
tween our techniques. HOD models describe galaxy bias us- 
ing a probability distribution P(N | M/j), the probability 
that a halo of mass Mh contains N galaxies of a given type. 
In HOD models, the satellite occupation function is con- 
structed using a parametric model. In contrast, in SHAM, 
the satellite occupation function is determined by the N- 
body simulation. Consequently, HOD models have several 
free parameters while the SHAM model used in this paper 
has none. 

In §2, we describe our methods - data, simulation, sub- 
halo abundance matching and rescaling technique. In §3, 
we describe our theoretical model and how it is sensitive 
to different cosmologies. In §4, we discuss our results and 
constraints on cosmological parameters. In §5, we discuss 
systematic uncertainties in our technique. In §6, we discuss 
the utility of our model in fitting other data sets to obtain 
complementary constraints on the cosmology. Finally, in §7, 
we summarise our results. 



2 METHODS 
2.1 Data 

Our galaxy clustering measurem ents are obtained fr om 7900 
deg^ of sky observed by SDSS (|Zehavi et all 1201 ll ). We re- 
strict our analysis to volume-limited samples of galaxies 

The clustering quantity we use for each sample is the 
projected auto correlation function Wp(rp) defined as: 



C(r-p,7r)d7r, 



(1) 



where rp is the projected separation between two galax- 
ies, TT is the line-of-sight separation between two galaxies, 
and (.{rpjir) is the measured two-dimensional correlation 
function. Due to the finite volume of the survey, the integral 
is limited to vTmax = 40 /i~^Mpc. 

Our analysis is primarily focussed on a volume-limited 
sample of galaxies brighter than Mr = -18.0. In addition, we 
make use of other volume-limited samples to test the robust- 
ness of our conclusions. The error estimates for each cluster- 
ing sample are com puted using the jackknife technique (see 
IZehavi et allboill for more details). 



2.2 Simulations 

We use two simu lations, the Millennium Simulation 
l|Springel et al.ll2005h and the higher resolution Millennium 



II simulation (jBovlan-Kolchin et al.|[2009l ) . Both simulations 
follow the evolution of 2160'^ particles from z — 127 to 
2 = in a ACDM cosmology (infiationary, cold dark mat- 
ter with a cosmological constant) with nM=0.25, f2A=0.75, 
h = J/o/100 km s~^Mpc~^= 0.73, primordial spectral in- 
dex ns = l, and the amplitude of mass fiuctuations , cr8=0.9 
where cts is the linear theory rms mass fiuctuation ampli- 
tude in spheres of radius 8 /i~^Mpc at 2 = 0. These pa- 
rameter value s were chosen to agree with WMAPl data 
ISpergel et all (|2003l ). and are different from but reason- 
ably close to current estimates from the cosmic microwave 
back ground (iLarson et ahl 1201 ll ) and large scale structure 
( Reid et al.ll2010l ). the main difference is that the more re- 
cent data favour a lower value of erg. 

The Millennium Simulation (MS) simulates a comov- 
ing box that is 500/i~^Mpc on each side while the Millen- 
nium II Simulation (MS-II) simulates a comoving box that 
is 100/i~^Mpc on each side. The simulation particle masses 
are rrip = 8.6x10* /i~^Mq in MS and 6.9x10^^ R-'^Mq in 
MS-II. 

For each output epoch in each simulation, the friends- 
of-friends (FOF) algorithm is used to identify groups by link- 
ing together particles sep arated by less tha n 0.2 of the mean 
inter parti cle separation (jPavis et al.lll"985h . The SUBFIND 



algorithm (jSpringel et al 



20011 ) is then apphed to each FOF 



group to split it into a set of self-bound sub halos. The cen- 
tral subhalo is defined as the most massive subunit of a 
FOF group. We construct subhalo merger trees which link 
each subhalo at each epoch to a unique descendent in the 
following epoch. These merger trees allow us to track the for- 
mation history of eac h (s ub)halo that is identifi e d at z — 0. 
ISpringel et all l|200lh and lBovlan-Kolchin et all (|200i) pro- 
vide a detailed description of these simulations and the post- 
processing techniques. 



2.3 Subhalo Abundance Matching 

Subhalo abundance matching (SHAM) is a technique for as- 
signing galaxies to simulated dark matter halos and subha- 
los. The essential assumptions are that all galaxies reside in 
identifiable dark matter substructures and that luminosity 
or stellar mass of a galaxy is monotonically related to the 
potential well depth of its host halo or subhalo. Some imple- 
mentations use the maximum of the circular velocity profile 
as the indicator of potential well depth, while others use 
halo or subhalo mass. The first cle ar formulations of SH AM 
as a systematic metho d appear in lConrov et al.l (|2006l ) and 
IVale fc Ostriked | 20061 ). but these build on a number of pre- 
vious studies that either test the underpinnings of SHAM 
or implicitly ass u me SHAM-like galaxy assignment (e.g. 
Colm et al.l I1999I : iKravtsov et"al] |2004| : iNagai fc Kravtsovl 



2001). 

N-body simulations produce subhalos that are located 
within the virial radius of halos. The present mass of sub- 
halos is a product of mass built up during the period 
when the halo evolves in isolation and tidal mass loss af- 
ter it enters the viri a l radius of a more mass ive halo (e.g. 
IKravtsov et al. I I2OO4I : iKazantzidis et al.ll2004l ). The stellar 
component, however, is at the bottom of the potential well 
and more tightly bound making it less likely to be affected 
by tidal forces. Therefore, several authors (e.g. lConrov et all 
l2006. : ,Vale fc Ostriker..2006 ) argue that the properties of the 
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stellar component should be more strongly correlated with 
the subhalo mass at the epoch of accretion rather than at 

z = Q. 

IVale fc Ostrikeil l|2006l ) apply a global statistical cor- 
rection to subhalo m asses r elative to halo mas ses (as do 
[Weinberg et al.ll2008l ). while IConrov et alJ (|2006t ) explicitly 
identify subhalos at the epoch of accretion and use the max- 
imum circular velocit y at that epoc h. Our formulation here 
is similar to that of IConrov et alJ (|2006), though we use 
mass rather than circular velocity. Specifically, we assume a 
monotonic relationship between galaxy luminosity and halo 
mass at the epoch of accretion and determine the form of 
this relation by solving the implicit equation 

ns{> Mr) =nH{> Mh), (2) 

where ns and tih are the number densities of galaxies and ha- 
los, respectively, is the galaxy r-band magnitude thresh- 
old, and Mh is the halo mass threshold chosen so that the 
number density of halos above it is equal to the number den- 
sity of galaxies in the sample. The quantity Mh is defined 
as follows: 

^ _ I A/halo = 0) for distinct halos, 
\M-tiaio{z — Zsii,t) for subhalos, 

where Zsat is the epoch when a halo first enters the virial 
radi us of a more mas sive halo. 

IGuo et"ai]||2012l ) compare the mass functions of SHAM 
selected subhalos (see their figure 1) in MS and MS-II finding 
good convergence over the mass range where both simula- 
tions have adequate resolution. At 2 = 0, they find that the 
(SHAM selected) subhalo mass function converges for sub- 
halos with more than ~ 10"^ particles at infall. This is more 
than an order of magnitude greater than the level at which 
the halo mass function converges (Bovlan-Kolc hin et al.l 
I2OO9I ) because subhalos experience tidal stripping after in- 
fall, and must therefore have a larger number of particles at 
infall to be reliably identified at z = 0. 

Fig. [1] shows the satellite occupation function as a func- 
tion of halo mass i.e. the mean number of SHAM selected 
satellite subhalos per parent halo as a function of parent halo 
mass. In both MS and MS-II, we rank order subhalos identi- 
fied at z = by Mh, mass at z = for central subhalos and 
at infall for satellite subhalos. We then select subhalos above 
a mass threshold of 8.6x 10^° /i"^M0 which is 100 times the 
particle mass in MS. Because of its higher resolution (factor 
of ~ 100 difference in particle mass), more satellite subhalos 
are identified in MS-II compared to MS across the range of 
parent halo masses. 

In MS, there is a missing population of satellite subha- 
los that have fewer than 10"^ particles at infall and cannot be 
identified at z — because of subsequent stripping, although 
similar objects are identified in the higher resolution MS-II. 
In order to exploit the larger box size and better statistics 
in MS, we construct a subhalo catalogue. Millennium Sim- 
ulation Plus (MS+), that matches the satellite occupation 
function of MS-II. To achieve this, we augment the subha- 
los identified at z = in MS with additional subhalos that 
are identified at high z, but cannot be identified at 2 = 0. 
These subhalos are randomly chosen from among disrupted 
subhalos that were above the mass threshold at infall. The 
number of additional subhalos in each parent halo mass bin 
is set so as to obtain a match to the MS-II satellite occupa- 




12 13 14 

log M„/h-'M<3 

Figure 1. Mean number of SHAM selected subhalos per halo 
as a function of parent halo mass at z = 0. Subhalos above an 
infall mass threshold of 8.6x10^'' h~^MQ are selected which cor- 
responds to 100 times the particle mass in the Millennium Simu- 
lation. 

tion function. The 2 = positions of each disrupted subhalo 
is obtained by tracking the position of its most bound par- 
ticle from Zsat to 2 = 0. 

The difference in the satellite occupation functions be- 
tween MS and MS-II depends on the mass threshold that is 
applied, and for a mass threshold above ~ IO^'^/i^^Mq, they 
converge. Therefore, the MS-I- catalogue has to be freshly 
constructed for each mass threshold that is required. While 
fitting the SDSS data, we carry out the procedure to con- 
struct MS-I- with a mass threshold which gives the galaxy 
number density of the observed sample. 

Fig. [2] shows the two-point correlation function of 
SHAM selected subhalos in MS-II compared to MS+. On 
scales not affected by box size effects, the auto correlation 
functions extracted from the MS-II simulation and our MS-I- 
method described above show strong agreement demonstrat- 
ing the efficacy of this procedure. 

2.4 Rescaling Technique 

Using the clustering of galaxies at 2 = to infer cosmologi- 
cal parameters requires modelling the growth of density per- 
turbations in the non-linear regime. In principle, one would 
need a numerical simulation for each set of cosmological pa- 
rameters. However, running such a suite of high resolution 
cosmological N-body simulations to adequately sample the 
cosmological parameter space is computationally expensive. 

Instead, we rescale our simulation from its "native" 
cosmology to different cosmologies with slightly different 
parameters, using a technique that involves rescaling the 
box length, velocity and mass units and relabelling output 
times. The first attempt to rescale the simul ation output 
to a d ifferent cosmology w as carried out by IZheng et al.l 
(|2002l) . lHarker' et al.l l|2007h used this rescaling technique in 
combination with semi-analytic models of galaxy formation 
to constrain the cosmological parameters. More recently, 
lAngulo fc Whitel l|2010l ) presented the technique of rescal- 
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Figure 2. Two-point correlation function of z = subhalos that 
are above an infall mass threshold. In the Millennium II simula- 
tion, we only include subhalos identified at 2 = 0, while in the 
Millennium simulation, we add a fraction of disrupted subhalos 
to those identified at z = so that we match the mean number 
of subhalos per parent halo as a function of parent halo mass in 
Millennium II (see Fig.[T]l. 



ing the output of a cosmological simulation from its "na- 
tive" cosmology to a cosmology with different parameters 
as a systematic algorithm. 

In this paper, we use the lAngulo fc Whitd (|2010l ) al- 
gorithm to generate (sub)halo catalogues for a given set of 
cosmological parameters. Following Ruiz ot al. (2011), we do 
not apply the final step of the lAngulo fc White! I|20ld ) algo- 
rithm which involves adjusting the amplitudes of large-scale 
linear modes which would in any case have negligible effects 
on the clustering statistics consider ed in this p a per. 

Angulo & White (2010) and iRuiz et all l|201ll ) find 
good agreement between the halo catalogues extracted from 
simulations scaled to a target c osmology a n d sim ulations 
run with the target cosmology. IGuo et all (|2012[ l extend 
this comparison to subhalos finding good agreement between 
subhalo mass functions extracted from rescaled and "native" 
simulations. 
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Figure 3. A pedagogical demonstration of the method employed 
in this paper. The effect on the two-point correlation function of 
SHAM selected subhalos of changing Qj^ at fixed erg in panel (a), 
and of changing erg at fixed f^M in panel (b). 



3 COSMOLOGICAL CONSTRAINTS FROM 
THE GALAXY TWO-POINT 
CORRELATION FUNCTION 

We probe the cosmology by fitting the observed galaxy two- 
point correlation function of a volume limited sample of 
galaxies. To generate a theoretical prediction for the two- 
point correlation function for a given cosmology, we imple- 
ment the following steps. First, we rescale our simulation to 
the target cosmology using the procedure described in §2.4. 
We then select (sub)halos using the SHAM procedure de- 
scribed in §2.3. Briefly, subhalos above a mass threshold at 
infall for satellite subhalos, and at z = for central subha- 
los are selected. Our subhalo mass threshold is set so that 
the number density of SHAM selected subhalos in our sim- 
ulation is equal to the number density of observed galaxies 
in the sample we are fitting to (see equation 2). Finally, 



we compute the projected two-point correlation function of 
these SHAM selected subhalos in the simulation cube. 

Fig.Oprovides a pedagogical demonstration of our tech- 
nique demonstrating that, in principle, it could be used to 
obtain cosmological constraints. Panel (a) shows the effect 
of changing SIm on the two-point galaxy correlation func- 
tion in a standard ACDM cosmology with all other cosmo- 
logical parameters held constant. Decreasing boosts the 
two-point correlat ion function on all scales (see figure 8 of 
IZheng et al. I I2OO2I for a similar result). Panel (b) shows the 
effect of changing erg on the two-point galaxy correlation 
function in a standard ACDM cosmology with all other cos- 
mological parameters held constant. Increasing erg boosts 
the two-point correlation function on all scales. 
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Figure 4. The solid curve is the galaxy two-point correlation 
function of our best-fit model with = 0.275 and erg = 0.86. The 
points with error bars are the SDSS observed galaxy two-point 
correlation function from a volume limited sample of galaxies with 
Mr < -18.0. 



4 COSMOLOGICAL CONSTRAINTS 

We determine constraints on the cosmological parameters by 
comparing the predicted two-point galaxy correlation func- 
tion for a grid of models in the ag,-O.M plane to the SDSS 
observed two-point galaxy correlation function of galaxies 
brighter than Mr = -18. We restrict our analysis to scales 
below 10 /i"^Mpc. 

We compute the for each model using the full co- 
variance matrix. The covariance matrix is calculated by 
(2011) using jackknife resampling (see §2.2 of 
[2OII1 for further details of the method). There 
are two free cosmological parameters, Q.m and erg, and no 
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other free parameters in our model. 

We impose the following flat priors on our cosmological 
parameters such that 0.2 < 57m < 0.35 and 0.65 < o-g < 1. 
Rescaling our simulation to a model with a different ampli- 
tude of clustering relies on relabelling the simulation output 
epochs. Since our simulations are only run to z = 0, we are 
restricted to erg < 0.9. For 0.6 < erg < 0.9, linearly inter- 
polating log(^) with erg works accurately, and so we extend 
this scaling to extrapolate the galaxy two-point correlation 
function to generate predictions for models with erg > 0.9. 

Fig. m shows the galaxy two-point correlation function 
for our best-fit model with = 0.275 and erg = 0.86 plot- 
ted against the data. The for our best-fit model is 9.21. 
As we have two free parameters and nine data points, this 
corresponds to a per degree of freedom of 1.3. The proba- 
bility of drawing a value > 9.21 from a distribution with 
7 degrees of freedom is 0.24, and so our best-fit model is an 
acceptable fit to the data. 

Fig. [5] shows our main result, cosmological constraints 
in the f^M-erg plane. Because of the way 57m and erg affect 
the two-point galaxy correlation function (see §3), our con- 
straints on them are correlated. Note that only the portion 
of the outer contour with erg > 0.9 depends on our extrap- 
olation of the two-point correlation function. 

Fig. |6] shows the marginalised probability distributions 




0.20 0.22 0.24 0.26 0.28 0.30 0.32 0.34 

Figure 5. Joint constraint in the trg - Qm plane. The inner 
contour shows the boundary of the 68% confidence region and 
the outer contour shows the 95% confidence region. 



of f^M and erg. Our constraints on the individual parameters 
are Q.m = 0.29 ±0.03 and erg = 0.86 ± 0.04 (68%). 

We emphasize that our constraints on these parameters 
are obtained for fixed values of other cosmological parame- 
ters. We defer discussion of the implications of this to the 
final section. 



5 SYSTEMATIC UNCERTAINTIES 

We discuss the systematic uncertainties in our analysis, fo- 
cusing on the uncer tainties in our theoretical modelling. See 
IZehavi et all lj201ll 'l for a discussion of observational uncer- 
tainties. 



5.1 Scatter in SHAM 

SHAM assumes a strictly monotonic relation between 
(sub) halo mass and galaxy r-band luminosity with zero 
scatter. While the assumption of zer o scatter is ide- 
alised, several observation al studies (e.g. IZheng et al. I l2007l : 
Ivan den Bosch et al.l l2007i ) indicate that the scatter in lu- 
minosity at fixed (sub) ha l o mas s is small. While from the 
theory side, ISimha et alj l|2012t) find a strong correlation 
between r-band luminosity and subhalo mass in their SPH 
simulations, with a small scatter of 0.15 dex in luminosity 
at fixed subhalo mass. 

Fig. [7] shows the effect of scatter in the (sub)halo mass 
- luminosity relation on the predicted gala xy two-point cor - 
relation function. We use the results of Si mha et al.l (|2012l ) 
to obtain an estimate of the scatter in subhalo mass in bins 
of galaxy luminosity, and perturb our subhalo masses by a 
random number drawn from a distribution centered on this 
value. We then compute the two-point correlation function 
for this sample of subhalos. We plot the fractional differ- 
ence between the two-point correlation function of this sam- 
ple and the fiducial model that assumes no scatter between 
(sub)halo mass and luminosity. For comparison, we show 
the fractional difference in correlation functions between our 
fiducial model and a ler change in erg calculated by taking 
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Figure 6. Marginalised posterior probability distribution of erg 
in panel (a), and of in panel (b) from fitting the two-point 
galaxy correlation function. 



the mean of the absolute value of the fractional difference 
in the correlation function from a positive and negative la 
change in erg from our best-fit value of o"8- The effect of ran- 
dom scatter in the subhalo mass - luminosity relation on the 
galaxy two-point correlation function is negligible compared 
to the effect of changing erg by one standard deviation. 

Fig. [7] shows the effect of scatter in the subhalo mass - 
luminosity relation for a particular set of cosmological pa- 
rameters and a particular realisation of the scatter. However, 
our investigations lead us to believe that the effect of scatter 
in the subhalo mass - luminosity relation is not sensitive to 
small changes in the cosmological parameters. 

Random scatter in the subhalo mass - luminosity rela- 
tion has negligible impact on the galaxy /subhalo two-point 
correlation function because the effect of such scatter is only 
seen at the boundary of our sample. Incorrect identification 
of a galaxy's host halo does not affect our result as long as 
the halo hosts a galaxy that is sufficiently luminous to be in- 
cluded in our sample. However, this effect could potentially 
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Figure 7. Effect of scatter in the subhalo mass - luminosity 
relation on the galaxy two-point correlation function compared 
to the effect of a Icr difference in a-g. Fractional difference in the 
two-point galaxy correlation function compared to the fiducial 
model is plotted against length scale. 



be a significant source of uncertainty for other statistics that 
depend on correctly identifying the host halos of galaxies. 



5.2 Satellite Galaxy Fraction 

SHAM relies on accurately identifying substructures and re- 
covering their properties. Since we trace the Zsat progenitors 
oi z = Q substructures, our results are unlikely to be af- 
fected by random fluctuations in the density field of halos 
that may be spuriously identifled as substructures. However, 
if subhalos hosting satellite galaxies that have merged with 
the central galaxy of the halo are identified as substructures, 
we would overestimate the halo occupation of massive halos. 
Conversely, subhalos that fall into more massive halos and 
lose a substantial fraction of their mass due to tidal strip- 
ping may no longer be resolved in the simulation at z = 
alth ough they m a y still host satellite galaxies. 

ISimha et al] (|2012D find that for a fixed subhalo mass, 
satellite galaxies in their SPH simulation are typically some- 
what less luminous than central galaxies because of differ- 
ences in the ages of their respective stellar populations. If 
this were to be true of the real Universe, SHAM would over- 
estimate the satellite galaxy fraction. 

Fig. [8] shows the effect of changing the satellite galaxy 
fraction. To increase the satellite galaxy fraction, we include 
additional "merged" subhalos i.e. subhalos that exist at high 
redshift, but have merged hy z = G using its most bound 
particle when it was last identified to track its position to 
z = Q. We then compute the two-point correlation function 
for this sample of subhalos. We plot the fractional differ- 
ence between the two-point correlation function of this sam- 
ple and the fiducial model. For comparison, we again show 
the fractional difference in correlation functions between our 
fiducial model and a Icr change in erg. 

Increasing the fraction of satellite galaxies increases the 
number of close pairs and thus boosts the two-point corre- 
lation function on small scales. A small increase is also seen 
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Figure 8. Effect of cfianging the fraction of satellite galaxies 
on the galaxy two-point correlation function compared to the ef- 
fect of a Icr difference in erg. Fractional difference in the two-point 
galaxy correlation function compared to the fiducial model is plot- 
ted against length scale. 

on large scales because highly biased objects in massive ha- 
les are present in a larger number of pairs compared to a 
sample with fewer satellite galaxies. 

An error of 10% in the satellite galaxy fraction would 
alter our best-fit estimate of ag by roughly half a standard 
deviation making this a more significant source of systematic 
uncertainty than scatter in the (sub)halo mass - luminosity 
relation. 



1.676, 1.12 and 0.318 10"^/i~^Mpc^ respectively. The sam- 
ple with Mr < -21.5 corresponds approximately to galaxies 
brighter than L* , the characteristic galaxy luminosity above 
which the number density of galaxies falls exponentially. 

For each of these samples, we generate the predicted 
galaxy two-point correlation function for our best-fit cosmo- 
logical model by measuring the clustering of SHAM selected 
subhalos above an infall mass threshold determined so that 
the number density of subhalos is equal to the number den- 
sity of galaxies in the sample. Each panel of Fig. [9] shows 
the galaxy two-point correlation function predicted by our 
best-fit model plotted against the observed galaxy two-point 
correlation function for the corresponding sample. Formally, 
our model is a good fit to the data for the volume limited 
samples with Mr > -18.5, -19 and -19.5, but not for the 
brightest sample with Mr > -20.5. 

The error bars in each panel show the diagonal errors 
on the observed correlation function. For our main sample 
and samples with relatively high number density, the error 
on the correlation function extracted from our simulation 
is negligible compared to the observational errors, and can 
therefore be ignored. However, because of our finite simula- 
tion volume, as we go to brighter samples with low number 
density, the errors on our predicted correlation function are 
sufficiently large to be comparable to the observational er- 
rors. 

When we fit models to these other galaxy samples by 
allowing the cosmological parameters to vary, we find that 
for each sample, the best-fit cosmological parameters are 
slightly different, but consistent with the constraints from 
our main sample {Air < -18). In principle, tighter constraints 
could be obtained from a joint fit to all the datasets, but this 
would require an estimate of their covariance. 



5.3 Rescaling Technique 

iRuiz et al.l dioiH ) compare the properties of halos in a 
rescaled simulation to a simulation run with the target cos- 
mological parameters. They find that over 99% of halos with 
more than 50 particles are recovered in the rescaled simu- 
lation. In the rescaled simulations, the masses of halos are 
systematically underestimated by ~5%. But because we use 
halo mass only to assign a rank to subhalos, we are not af- 
fected by this systemic bias. However, it would be a source 
of error if a statistic that depended directly on halo mass 
were being used. The rescaled halo mass - "native" halo 
mass relation displays scatter. For our purpose, it would 
mimic the effect of scatter in the subhalo mass - luminos- 
ity relation discussed earlier. The positions of halos in their 
rescaled simulations are recovered to a precision of greater 
100 fe-^kpc. 

5.4 Other Galaxy Samples 

Our cosmological constraints have been obtained by fitting 
our model to the observed clustering of one volume-limited 
sample of galaxies with Mr < -18. As a consistency check, 
we compare the clustering predictions of our model with 
four ot her volume limite d samples of galaxies taken from 
SDSS ijZehavi et al.ll201ll l with Mr < -18.5, -19, -19.5 and 
-20.5 corresponding to a number density of galaxies of 2.311, 



6 COMPLEMENTARY CONSTRAINTS ON 
COSMOLOGY 

As we discussed earlier, (see Fig. [SJ , our constraints on Qm 
and £78 are correlated because of the way the changes in these 
parameters affect the galaxy auto-correlation function (see 
Fig. [31). To a certain extent, the effect of a higher Qm can 
be compensated for by a higher ag if only the galaxy two- 
point correlation function is used to constrain the cosmology. 
But complementary constraints on the cosmology can be 
obtained by probing other observables for which predictions 
can be generated using the same SHAM model used in this 
paper. 

Fig. [To] shows the effect of the cosmological parameters 
on the mean number of satellite galaxies per halo, {N}^^^^ as 
a function of halo mass in our simulation using our rescaling 
technique and SHAM model. The number of satellite galax- 
ies in a given halo, A'aat = A'^ - 1 where A'' is the number of 
galaxies in the halo. Panel (a) shows the effect of changing 
Qm at fixed as- Changing SIm affects the halo mass function 
making halos with a given number of satellites more massive 
for higher Qm- Panel (b) shows the effect of changing erg at 
fixed Qm- Increasing erg increases the number of halos of a 
given mass. Therefore, for a fixed number density of galax- 
ies, there must be fewer galaxies in each halo compared to 
models with lower erg. 

Decreasing erg and decreasing Qm, both boost {N)^^^ 
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Figure 9. In each panel, the points with error bars are the 
SDSS observed galaxy two-point correlation function in a volume 
limited sample of galaxies brighter than = -18.5, -19, -19.5 
and -20.5. The solid curve in each panel is the galaxy two-point 
correlation function predicted by our best-fit model with S^m = 
0.275 and erg = 0.86 for the corresponding galaxy sample. 



rather than counteracting each other as they do for the 
two-point correlation function. Moving along the degeneracy 
curve in the JIm-cts plane in Fig. [5] would generate different 
and easily distinguishable distributions of {N)^^_^. Therefore, 
for a given cosmology, the additional requirement of match- 
ing the galax;y two-point correlation function strongly con- 
strains the HOD, and consequently a direct measure of the 
HOD would place constraints on the cosmology that are 
complementary to our constraints from fitting the two-point 
correlation function. 

One possible direct measure of the HOD is the mean 
number of galaxies in halos of a gi ven mass. The ratio of 
these quantities, M/N, is used by [Tinker et all (|2012l ') to 
place constraints on cosmological parameters. 



7 DISCUSSION AND CONCLUSIONS 

We have placed constraints on erg and flu by comparing 
the SDSS observed projected galaxy two-point correlation 
function for a volume-limited sample of galaxies with Air < 
-18 to our model predictions generated using N-body simula- 
tions rescaled to the tar get cosmology using the technique of 
lAngulo fc Whitel l|2010l ) and populated with galaxies using 
subhalo abundance matching (SHAM). 

Assuming a fiat ACDM cosmology with ns~l, we find 
f^M = 0.29 ±0.03 and erg = 0.86 ± 0.04 at 68% confidence. 

Fig. [TT] compares our c onstraint in the flm- o's plane to 
constraints from WMAP7 jKomatsu et al.ll201ll '). Our esti- 
mates of both as and f^M are high compared to WMAP7, 
but are consistent at the ~2 a level. 

Our constraints are obtained for fixed values of other 
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Figure 10. Panel (a) shows the effect of changing Qm at fixed 
(78, and panel (b) shows the effect of changing erg at fixed on 
the mean number of satellite galaxies per halo as a function of 
halo mass. 



cosmological parameters. In contrast, the WMAP7 con- 
straints are obtained by marginalizing over all other cos- 
mological parameters. One significant implication of this is 
in regard to the shape of the power spectrum. We have as- 
sumed that the value of the primordial spectral index, ns = l 
in our simulation. In contrast, the WMAP-7 best-fit value 
of ns is 0.96, and ns=l is excluded at more than 2 a. While 
we are unable to comment on the effect of setting ns = 0.96 
on our constraints, forcing ng = 1 while fitting the CMB 
data would result in a higher best-fit value of erg. 

Our results a re consistent with and comparable to 
[Tinker et all l|2012l ) who fit the SDSS galaxy two-point cor- 
relation function and M/N (cluster mass to number ratio) 
using their HOD models, finding ag = 0.85 ± 0.06 and Q.m 
= 0.29 ±0.03. 

Ho wever there is s ome t ension between our results and 
those of iHarker et al] (|2007D who use semi-analytic models 
to populate N-body simulations rescaled to a given cosmol- 
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Figure 11. Joint constraint in the erg - Cm plane. The inner 
contour shows the boundary of the 68% confidence region and 
the outer contour shows the 95% confidence region. The filled 
contour is the result from t his work while the b lack solid open 
contours are from WMAP7 llKomatsu et al.ll201lll . 



ogy using a technique similar to that of this paper, and fit 
to the SPSS clustering find ing as = 0.97 ± 0.06. In con- 
trast to iHarker et all (j2007l ) who use a semi-analytic model 
of galaxy formation to populate their N-body simulation 
with galaxies, we use subhalo abundance matching which 
only assumes a monotonic relationship between galaxy lu- 
minosity and subhalo mass at infall. Secondly, the resolution 
of MS-II is ~ 2000 tim es higher than the sim u lation s used 
bv lHarker et all (j2007l ). Thirdly. iHarker et al.1 (|2007l ') use a 
Monte Carlo scheme to generate a merger tree for a halo 
based on its mass. Consequently, satelhte galaxy positions 
are not obtained from the simulation, and they are instead 
placed on random particles within the halo. We are unable to 
quantify the effects of e ach of t hese f actors. However, repeat- 
ing the work of Harker et al.l ()2007h with a high resolution 
N-body simulation and merger trees and subhalo positions 
extracted from the N-body simulation would be interesting 
and could potentially reveal the source of the tension. 

Besides the CMB and galaxy clustering data, several 
other methods have been employe d to constrain erg. For ex- 
ample, iMandelbaum et al.1 (|2012l ') find ffg(nM/0.25)° '^^ = 
0.80 ± 0.05 using galaxy-galaxy lensing, ISehak et al.1 ()2005l l 
find (Tg = 0.90 ± 0.03 by combining their analysis of the 
Lyman a forest power spectrum with CMB results. 

Because our simulation is only run to 2; = 0, we are 
unable to rescale our subhalo catalogue to cosmologies with 
erg > 0.9. Our model predictions for erg > 0.9 are based on 
extrapolating the correlation function as a function of the 
cosmological parameters. While this appears to be a rea- 
sonable approximation, it is likely that it has shortcomings 
that will be exposed by a simulation that is run with a higher 
value of (Tg. 

We have examined some potential sources of systematic 
error in §3 finding that random scatter in the luminosity - 
halo mass relation does not affect our results significantly. 
However, if the scatter were to be correlated with other prop- 
erties of the galaxy population, our results would be affected 
differently. We also examine the impact of a systematic er- 



ror in the fraction of satellite galaxies which could either 
be underestimated or overestimated in the models either 
for numerical reasons relating to identification of subhalos 
or due to systematic differences in the growth of central 
and satellite galaxies that violate the implicit assumptions 
of SHAM. We find that a 10% change in the fraction of 
satellite galaxies would alter our constraint on as by half a 
standard deviation. The two-point galaxy correlation func- 
tion on small scales is sensitive to the distribution of galaxies 
within massive halos. Therefore, our results would be signifi- 
cantly affected if the distribution of galaxies in halos were to 
systematically differ from the distribution of SHAM selected 
subhalos in N-body simulations. 

Despite these caveats, we emphasize that the technique 
presented in this paper can provide tight constraints on the 
cosmology using only low z data. The remarkable tightness 
of our constraint arises from the fact that unlike statistical 
descriptions of the distribution of galaxies in halos provided 
by HOD models or the CLF (Conditional Luminosity Func- 
tion), we do not have the freedom to define the HOD. Fur- 
thermore, our model does not make any assumptions about 
galaxy bias or the detailed physics of galaxy formation ex- 
cept for requiring a monotonic relationship between galaxy 
luminosity and subhalo mass at infall. Additionally, we as- 
sume that galaxies are the observational counterparts of sub- 
halos identified at 2; = in an N-body simulation, and that 
each observed galaxy can be associated with a subhalo. 

Although the potential for statistical tightening of the 
parameter constraints is limited, substantial improvements 
in the robustness of this technique can be achieved by future 
investigations. Firstly by adopting the shape of the power 
spectrum P{k) inferred by the latest CMB observations. Sec- 
ondly, a large volume simulation with the resolution of Mil 
or higher would remove the need for the correction to the 
subhalo fraction that we apply in this work and the uncer- 
tainties associated with it. Systematic uncertainties associ- 
ated with the rescaling technique scale with the magnitude 
of the difference in cosmological parameters between the 
"native" and rescaled cosmology. To minimise this, it would 
be useful to carry out a suite of simulations with different 
cosmological parameters that span the requisite range, and 
only apply the rescaling technique to generate predictions 
for intermediate values of the parameters. Finally, further 
investigations of the SHAM technique, particularly with re- 
spect to scatter in the luminosity - halo mass relation, and 
the distribution of subhalos within halos will help to clarify 
potential sources of systematic errors involved in this tech- 
nique. 
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